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Abstract 

We study the final-state interaction (FSI) effects in charmless B^^d^s PP decays. We consider 
a FSI approach with both short- and long-distance contributions, where the former are from in- 
elastic channels and are contained in factorization amplitudes, while the latter are from the residual 
rescattering among PP states. Flavor SU(3) symmetry is used to constrain the residual rescattering 
5-matrix. We fit to all available data on the CP-averaged decay rates and CP asymmetries, and 
make predictions on unmeasured ones. We investigate the Kir direct CP violations that lead to the 
so-called Ktt puzzle in CP violation. Our main results are as follows: (i) Results are in agreement 
with data in the presence of FSI. (ii) For B decays, the tt^tt" and Tr'^vr" rates are suppressed and 
enhanced respectively by FSI. (iii) The FSI has a large impact on direct CP asymmetries („4) of 
many modes, (iv) The deviation(A^) between A{B^ K^ir^) and A{B~ K~tt^) can be 
understood in the FSI approach. Since A{K~7r^) is more sensitive to the residual rescattering, the 
degeneracy of these two direct CP violations can be successfully lifted, (v) Sizable and complex 
color-suppressed tree amplitudes, which are crucial for the large tt^tt^ rate and A^, are generated 
through exchange rescattering. The correlation of the ratio B{-K^7r^)/ 8(71^71^) and is studied, 
(vi) The B~ 7r~-K^ direct CP violation is very small and is not affected by FSI. (vii) Several 
Bs decay rates are enhanced. In particular, the r]'r]' branching ratio is enhanced to the level of 
1.0 X 10^^, which can be checked experimentally, (viii) Time-dependent CP asymmetries S in B^^g 
decays are studied. The AS{B^ Ksf]') is very small (< 1%). This asymmetry remains to be 
one of the cleanest measurements to search for new physics phases. The asymmetry S from Bs to 
PP states with strangeness S= +1 are expected to be small. We found that the IS"! for B^ rjrj, 
ijrj' and rj'ri' decays are all below 0.06. CP asymmetries in these modes will be useful to test the 
SM. 

PACS numbers: ll.30.Hv, 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 



The study of B decays provides many useful information of the flavor sector of the 
Standard Model (SM) In particular, the measurements of the time-dependent CP 

asymmetries in kaon and charmonium final states give a rather precise value of sin 2/3 = 
0.681 ±0.025 y, where (3/(f)i = a-TgiV;^) with V the Cabbibo-Kobayashi-Mashikawa (CKM) 
matrix. In the SM, time-dependent CP asymmetries in penguin dominated modes are ex- 
pected to be close to the sin 2/5 value [SI]. Since the penguin loop amplitudes are sensitive 
to high virtuality, new physics beyond the SM may contribute to the time dependent CP 
asymmetries through the heavy particles in the loops. Consequently, these asymmetries are 
promising places to search for new physics effects s], 0, S, S Q, Is]. 

The measurements of direct CP violation {A) in B decays are also very useful and in- 
teresting. The A{B^ K^TT^) asymmetry was the first measured direct CP violation 
in B decays. The data confirmed a large A{B^ K~7t~^) with a negative sign as pre- 
dicted in perturbative QCD (pQCD) calculations |9|. On the contrary, although A(B~ 
K^n^) ^ A{B'^ K^TT^) was expected in many early theoretical predictions [ol, 3], the 
experimental evidence has been accumulated favoring a positive A{B~ —>■ K'n^"^ The 
recent measurements show A{K^7i^) = 
ing AA{Kn) = AiR-n^) - A{K-'k+] 



-9.8 



and A{K-7i^) = (5.0 ± 2.5)% [2], giv- 
;i4.8t|^)%, which is more than 5 cr from zero. 
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e in direct CP violation, which has attracted a lot of atten- 



19| . Several suggestions were put forward to resolve this 



puzzle. For example, some authors introduced next-to-leading order contributions in factor- 
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ization amplitudes [iJ], while some suggested new physics origins 
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for the deviation. 

It is well known that we need both weak and strong phase differences to have a non- 
vanishing direct CP violation. Strictly speaking the final state interaction (FSI) is the only 
source for non-vanishing strong phases. In addition, it is capable of enhancing the decay 
rates of many modes, which are measured to be larger than expected. For example, the 
large observed B^ ti^tt^ rate, which remains puzzling and is still posing tension in many 
theoretical calculations, can be obtained by using FSI 20|. Furthermore, it was realized 
recently that long-distance FSI may play an indispensable role in charmful as well as in 
charmless B decays 
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Data for Bg decays are starting to emerge from the Tevetron \^ and from B factories, 
and we anticipate more to come in the near future, from LHCb and other LHC experiments. 
Measurements of rates and CP asymmetries in Bg decays will be useful in testing the SM 
and in searching for new (physics) phases. In fact, recently, a claim on the evidence of new 



physics effect in the Bg mixing was put forward [23j. It is thus timely to study Bg decays. 

In this work, we investigate the effects of FSI on all charmless Bu,d,g —>■ PP decay rates 
and CP asymmetries. We outline the underlying physical picture of the FSI approach 
employed in this work. The master formula of FSI for charmless B —>■ PP decays is (see 
appendix A, if a derivation is needed) 

N 

Ar=j:^ie^. (1) 

k=l 

where A^^^ and are B decay amplitudes with and without FSI ^, respectively, i = 1, . . . ,n, 
denotes all charmless PP states, k = 1, . . . ,n,n + 1, . . . , N, denotes all possible states 
that can rescatter into the charmless PP states and S is the strong interacting S-matrix. 
Note that no approximation has been made in the above equation, which, in principle, all 
charmless B decay amplitudes should follow. In practice, this master formula is hard to use 
as it involves many states (the number N is in general quite large in a typical charmless B 
decay) . 

Let us investigate further the difficulties of using the above master formula. The number 
of states allowed to enter the formula grows with the mass of the decaying particle. For 
a typical B decay, there is a large number of the states involved in the equation and the 
contributions are hard to handle. For example, we may need to consider a rescattering 
process contributed from a multi-body final state, where the decay amplitude and the cor- 
responding rescattering ^-matrix element are poorly known. Therefore, the complication 
originates from the largeness of ms- However, it is precisely the largeness of rri R th at makes 
factorization approaches such as pQCD j^, QCD factorization (QCDF) [lO, 24] and soft 



collinear effective theory (SCET) [25| possible. These approaches achieve accessibility and 
simplifications. The underlying reason for the simplicity is related to the so-called duality 
argument, which uses the fact that when contributions from all hadronic states at a large 
enough energy scale are summed over, one should be able to understand the physics in 



^ Note that A^^^ contains weak as well as strong phases, while A'^ only has weak phases. 
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terms of the quark and gluon degrees of freedom. Hence, it is reasonable to expect that the 
main effect of FSI, especially those from inelastic channels, is included in the factorization 
amplitudes - a statement we expect to hold perfectly in the rrif, oo case. Since in the 
real world m^, is finite, whether it is large enough to validate the above argument should be 
answered by experiments. 

It is fair to say that most factorization results on CP-averaged charmless B PP decay 
rates, especially color-allowed ones, agree well with data. However, some measurements seem 
to imply the needs of sub-leading contributions. For example, rates of some suppressed decay 
modes, such as the above mentioned 5° n^n^ rate, and some CP-odd quantities, such 
as the K^tt^ direct CP violation, do not agree well with predictions. These are 

places, where sub-leading effects, such as FSI, could be visible. Therefore, although we 
expect factorization amplitudes to contain most of the FSI effects demanded in Eq. ([1]), it 
is likely that residual rescattering is still allowed and needed in B ^ PP decays at the 
physical tub energy scale. The group of charmless PP states is unique to the processes we 
are studying and is well separated from all other states. Since the duality argument cannot 
be applied to these states of limited number, part of their FSI effects may not be included 



1,0- 



in the factorization amplitudes |26| 

In summary, FSI in B decays may be simpler than we thought, since uib could be large 
enough to apply a factorization approach for the main part of FSI contributions. We may 
only need to include the left-over FSI, i.e. residual rescattering, in addition to the short- 
distance FSI in the PP sector. In this sense, the FSI approach we are using is a mild 
extension to the factorization approaches. 

Note that a similar approach analyzing early data was used in 27|]. There is one major 



difference. In 



271], in principle, no short-distance phase was allowed in factorization am- 



plitudes to avoid double counting, while here we do need short-distance phases to account 
for the FSI effects from all in-elastic plus some quasi-elastic channels. There are also other 
works in the literature discussing rescattering among PP states and/or from some in-elastic 



channels 
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29, 3^, 
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32|]. For example, in 22], rescattering from PP and 



final states was considered, and the main FSI contributions resemble the charming penguin 
ones 
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331 ]. We also note that similar discussion of the factorization of S into short-distance 



and residual parts, as well as the discussion of the approximation done when assuming 5* is 
block diagonal (with a block for the PP states) can be found in \2^ . 
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The layout of the present paper is as follows: In Sec. II we introduce the formalism. We 
then use it in Sec. Ill to study Bu,d,s PP decays. Results and discussions are presented in 
Sec. IV. Sec. V contains our conclusions. Some derivations, including those lead to Eq. ([T]), 
are given in Appendices. 



II. FORMALISM 



In this section, we develop the formalism. The reader who is not interested in the detail 
of the formalism may proceed directly to the numerical analysis section. 
Without loss of generality, we can re-express the S'-matrix in Eq. ([T]) as 

n 

Sik = ^^('5i)ij(iS2)jfc5 (2) 

i=i 

where Si is a non-singular x n matrix with n the total number of charmless PP states and 
^2 is defined through the above equation, i.e. ^2 = S^^S. The physical picture mentioned 
in the last section is close to the one in factorization approaches, except that there are still 
some residual rescattering effects, and we have 

N 

Si = Sres, ^j"^^ = ^^('^2^ (3) 

k=l 

with the total number of states entering Eq. ([1]) and the factorization amplitude. The 
residual rescattering effect is encoded in the Sres niatrix. Note that although S is unitary, 
Sres needs not be so, as it describes the residual rescattering among various charmless PP 
states. In factorization approaches, the above Sres is taken to be unity. We shall use the 
up-to-date data to determine Sres- It should be reminded that our framework does not 
exclude the fully factorized case {Sres = 1) and, hence, it is also being tested. To apply the 
above formula, we need to specify the factorization amplitudes. In this work, we use the 
factorization amplitudes obtained in the QCD factorization approach 24 1. 
Combining Eqs. ([1]) and ([3]), we have 

n 

= J2isu^hAf-, (4) 

where, as mentioned before, i,j = 1, . . . ,n denote all charmless PP states. The number of 
parameters needed to describe Sres seems to be quite large. This is, however, not the case. 



since strong interaction has (an approximate) SU(3) symmetry, which is expected to be a 
good one at the ttib rescattering scale and, hence, can be used to constrain the form of Sres- 
Exphcitly, through SU(3) symmetry, we have 
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where we have i5^ef j = (1 + iTj)^/"^, with 
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and 



= diag(ro, ro, tq, ro, fo, ro, ro, ro, ro) 
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(10) 

The rescattering parameters ro,a,e,t, ro,a,e,t, ro,a,e,t, ro,a,e,t and fo,a,e,t denote rescattering in 
n(8)n(8) n(8)n(8), n(8)n(8) ^ n(8)r7i, n(8)r/i ^ n(8)r/i and r/ir/i r^ir^i, re- 
spectively, and the subscripts 0,a, e,t represent flavor singlet, annihilation, exchange and 
total-annihilation rescatterings, respectively (see Fig. [1]). Note that for identical particle 
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FIG. 1: Pictorial representation of (a) charge exchange r^, (b) singlet exchange tq, (c) annihilation 
and (d) total-annihilation rt for PP (re)scattering. 

final states, such as tt'^tt", factors of are included in the amplitudes and the corre- 

sponding Sres matrix elements. The Prjs, Prji are not physical final states. The physical 
7], 7]' mesons are defined through 







[r/l 





COS — sin 
sin cos ^ 

with the mixing angle 'd ~ —15.4° js^. For the r/(')r/(') states, we have 

^ cos^-i? — -\/2 cost? sin -(9 sin^ ^ ^ 

a/2 cos -(9 sin "i? cos^ -i? — sin^ -(9 — v^cos'i^sin-i? 
sin^ "(9 cos sin cos^ 



(11) 



riT]' 



\VV ) 



\ 



f 



VsVs 
VsVi 
\ViVi j 



where the identical particle factor of 1/ \pl is properly included in the mixing matrix. 

The matrices 71^2,3,4 can be obtained through a diagrammatic method by matching 
the Clebsh-Gordan coefficients of scattering mesons (see Fig. [1]) or by using an operator 
method. We have TriJX^^l^W^^X^^^W^^) jl, Tr (ni'^n'2°nf*n°"*), Tr(ni°n°"')rr(n!fn°"*) and 
Tr(n™n2°)Tr(n°"*n2"*) corresponding to rg, r^, tq and rt contributions, respectively (see 
similar discussion for the case of charmful B decays in 2l|). Note that due to Bose- Einstein 
statistic and the S-wave configuration of the final state mesons in i? ^ -P1-P2 decays, the 
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amplitude should be symmetric under the exchange of the indices 1 and 2. Consequently, the 
above terms exhaust all possible combinations for 11(8) n(8) — > n(8) n(8) scatterings. For 
operators involving rji, we suitably replace 11 in the above expressions by r^ilsxs to obtain 
operators corresponding to fj, fj, fj and f j. 

It can be easily seen that rescattering formulas for charmless Bg PP decays resemble 
those for — > PP decays. Information on Srec obtained from 5° decays can be used to 
predict Bs decay rates. 

At first sight, it appears that we need 40 real parameters (from 20 complex rescattering 
parameters: ro,a,e,i, ^o,a,e,t, ^o,a,e,t, ^o,a,e,t and fo,a,e,t) to describe Sres- The number of inde- 
pendent parameters is actually much lower for two reasons. First, rescattering parameters 
enter Sres only through 7 independent combinations: 1 + Oj-q + r^), z(re — r^), iija + rt), 
i{2fa + fe), 1 + + {Afa + 2fe)/3], i(4f« + 2f^ + 3ft) and 1 + i[fo + {Afa + 2fe + 3n)/6]. 
Second, SU(3) symmetry imposes further constraints on these combinations. 

Flavor symmetry requires that {Sres)^ with an arbitrary power of m should also have the 
same form as Sres- More explicitly, from SU(3) symmetry, we should have 



{SresT = {l + lTr = l + lT^'^\ 

where T*^"^) is defined through the above equation and its form is given by 
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0, a, e, t. 



It is found that the solutions to Eqs. (fT3!) and (fT^ are given by 
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where UJ^ and VJ" are elements of 
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— sm u cos ^ / \ U e i / \ sm u cos u 

respectively. From the above solution, we see that two real mixing angles r and z/, and five 
comp/ex phases 527,8, ^'s i needed to describe (Sj-es)"^ in the full SU(3) case. 

Several remarks are in order, (i) The subscripts of phases denote the corresponding SU(3) 
multiplets and more details will be given shortly, (ii) The imaginary parts of 527,8,i7<^8i 
control the leakage of the non-unitary through the scattering of PP states into non-PP 
states, (iii) As we shall see, the iS™^ can be factorized into two parts depending only on the 
real and the imaginary parts of these phases, respectively, (vi) To reduce the number of the 
FSI parameters we will consider a restricted SU(3) case, which is close to a U(3) symmetric 
case. 

Since charmless mesons P consist of an octet (n(8)) and an singlet (771), we have 8 8, 
8 ® 1, 1 ® 8 and 1 ® 1 SU(3) products for P1P2 final states. Due to the S'-wave configuration 
of P1P2 in B decays and the Bose-Einstein statistics, the resulting SU(3) multiplets should 
be symmetric under the exchange of Pi and P2. The allowed ones are the 27, 8 and the 1 
Tom 8 ® 8, the 8' from the symmetrized 8 ® 1+1 8, and 1' from 1 ® 1 (see, for example 



351]). Hence, from SU(3) symmetry and the Bose-Einstein statistics, we have 



27 



(5,e.r = 5^|27;a)e^™'^-(27;a| + ^ \p;b)U^^{q;b\ + b; l)Vp^(g; 1|, (17) 

a=l b=l p,q=8,8' P,g=l,l' 

where a and b are labels of states within multiplets. It can be easily seen that the above 
form of is preserved for any value of m. We note that similar formulas for B PP 



rescattering (excluding P = rji) from SU(3) symmetry have been used in 271. |28[| 



From Eq. (JT71) we see that the matrix is in block-diagonal form and we also have 

V"^(z/,5i,5;) = V'"(z/,Re(5i,Re5;) ■ V'"(z/,zlm(5i,zlm(5;), (18) 

which can be proved by using the explicit expressions of and V"* given in Eq. (fT6ll . 
Consequently, the matrix can be factorized into two matrices containing only real and 
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imaginary phases, respectively, i.e. 

'^1,8,27, 5'i,8) = '5™,(r, i^;Re5i,8,27, Re^i s) ■ 5™,(r, z/; zlm5i,8,27,iilm5i 8). (19) 
Note that iS^s(t, z/; i Im i Im (5-) is a n x n real matrix. Substituting the above expression 

1 /2 

of S'res into Eq. (jl]), we have 

^F57 = Re 5i,8,27, Re 5i g) ■ S^^^{t, u; i Im ^1,8,27, i Im §) " ^fac- (20) 

An overall phase in Eq. fl20l) can be removed and we are free to set Re ^27 = 0. Furthermore, 
in our analysis (as well as in many analyses using naive or QCD factorization approaches), 
various form factors and rris in A/ac are allowed to float in some given ranges of values. 
Therefore, an overall scaling factor (exp(— Im527)) can be absorbed into the form factors in 
Afac and we set Im527 = to avoid double counting. We are left with two mixing angles, 
four real phase differences and four imaginary phase differences: 

T, z/, = ReiSsio - S27), o^'^ = Re((5i(,) - ^27), 

= Im(58„) - ^27), i^'^ = Im(5i(,) - ^27). (21) 

The number of the residual FSI parameters is still quite large. It will be preferable to reduce 
it through some physical arguments or the consideration of some plausible cases. 

It is interesting to see how the residual FSI behaves in a U(3) symmetric case. It is known 
that the 11^(1) breaking is responsible for the mass difference between 77 and rj' and U(3) 
symmetry is not a good symmetry for low-lying pseudoscalars. However, U(3) symmetry 
may still be a reasonable one for a system that rescatters at energies of order uib- The mass 
difference between rj and r/', as an indicator of U(3) symmetry breaking effect, does not lead 
to sizable energy difference of these particles in charmless B decays. In the literature, some 



authors also make use of U(3) symmetry in charmless B decays (see, for example 36l|). 
The full U(3) symmetry requires: 

ri = fi = fi = fi, (22) 

for each i = 0,a, e, t. This imposes a major reduction of parameters. Note that the reduction 
is more easier to preform in the r, formalism than in the SU(3) decomposition formalism. 
This is one of the advantages of the former formalism. 
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In the U(3) case, we are constrained to have (see Appendix B) 

^M^M _ Q_ ^23) 

Consequently, there are two different solutions: (a) the annihilation type (ri™'' 7^ 0, ri"^^ = 0) 
with 

S27 = 6'^ = 6[, 68, 61, T = 
and (b) the exchange type (re™"* 7^ 0, ri"^^ 

^27 = ^8 = '^1' '^8 = ^1, T 

The explicit expressions of r^"^^ in terms of these phases can be found in Appendix B. 
It is interesting to note that in both solutions of the U(3) case, a common constraint 

527 = ^8 = S[, (26) 

has to be satisfied. To reduce the number of the residual FSI parameters shown in Eq. (l2Til . 
we consider a restricted SU(3) case, which is close, but not necessarily identical, to the full 
U(3) case. Motivated by Eq. (|26|) . we consider the parameter space around 

5' ~ a ~ 0, ~ ^' ~ 0. (27) 

The above restriction on the FSI parameter space is a rather strong model assumption. 
When comparing the fitted FSI parameters with those in Eqs. (IMI) and (1251) . it is possible to 
determine whether the exchange-type, the annihilation-type or a mixed solution is preferred 
by data. 



1 . _i4V5 

- - sm , 

2 9 ' 



1 . _i4v^ 
--sm , 

2 9 ' 



(24) 



^(m) 



0) with 



1 . _i 

- sm , 

2 3 ' 



1 „i 2^2 
z/ = - sm . 

2 3 



(25) 



III. NUMERICAL RESULTS 



In our numerical study, masses and lifetimes are taken from the review of the Particle 

n 

Data Group (PPG) p|, and the branching ratios of B to charmless meson decays are taken 



from 



0.0 



38|. We use = 131 MeV, Jk = 156 MeV [l] and /s^^, = 200 (230) MeV for 



decay constants. The values of CKM matrix elements are taken from the central values of 
the latest CKM fitter's results [39 1. 
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TABLE I: Confidence level (C.L.), Xmin/'^-^-^- ^^'^ various contributions to Xmm best fitted 



solution. Numbers of data used 


are shown in 


parentheses. 






C. L. Xmi„./d-o.f. 




.} ^{A(B«^Kt!), 


..} X{8(B-_^_ft:7r),... } 


A{^(B-^K7r),...} 


0.04 (43) 1.51 (43) 


4.6 (4) 


1.6 (3) 


4.5 (4) 


7.0 (4) 


2 2 
^{e(B--»7r7r),... } X{A{B-^7nr),...} 


^{B(B0^7r7r),. 


. } ^{A(B°->TnT), 




As(B«))} 


2.6 (4) 6.8 (4) 


8.0 (9) 


2.4 (3) 


1.6 (4) 


6.0 (4) 



We use the QCD factorization calculated amplitudes 2J] for the factorization amplitudes 
in the right-hand-side of Eq. (jl]). We take the renormalization scale /i = 4.2 GeV and the 
power correction parameters Xa,h = ln(mB/A/i)(l + pa.hg^'^^'")- Hadronic parameters in 
factorization amplitudes are fit parameters in addition to FSI parameters, and are allowed 
to vary in the following ranges: 

< Pa = P// < 2, -TT < (t)A,H < vr, m3(2.1GeV) = (100 ± 30)MeV, 
F(f^(0) = 0.25 ± 0.05, F(f^(0) = 0.35 ±0.08, F^^^(O) = 0.31 ± 0.08. (28) 



Note that we take pa = Ph for simplicity. These estimations agree with those in [2J, |4l|, |42| , 
while the ranges of form factors are slightly enlarged to include the possible effect of the 

1 /2 

overall scaling factor exp(— Im527) from SrL ■ Other parameters (if not specified explicitly) 



in the QCDF amplitudes are taken from the central values of those used in 2J|. For the 
FSI parameters, we set allowed ranges to be 

-|<r,z/<^, -TT < (5, a < TT, -0.35 < K,^ < 0.35 (29) 

for the mixing angles, real and imaginary parts of FSI phase differences. In the fit we take 
6' = a' = k' = C,' = as mentioned in the end of the previous section. The effects of relaxing 
these constraints will also be estimated. 

We perform a analysis with all available data on CP-averaged rates and CP asymme- 
tries in Bu^d,s PP decays. There are altogether 43 data used in the fit. The confidence 
level and for the best fitted case are shown in Table [B Contributions to Xmin froiri vari- 
ous sub-sets of data are also given. For example, X|g(;BO^X7r) } table denotes the 
contribution obtained from 4 CP-averaged — > K^tx^ , K^tt'^, K^t], K^t]' rates, which are 
related through FSI [see Eq. ([5]), and see Eqs. ([6])-([8]) for other groups]. Numbers of data 
used are shown in parentheses. 
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TABLE II: Fitted hadronic and FSI parameters. Upper table contains fitted parameters in fac- 
torization amplitudes, while the lower one contains fitted FSI parameters. Note that parameters 
with values given in parentheses are not fitted ones (see text). 



PA,H 






m^(MeV) 


Fo^-(O) 








-■-■-^0-0.23 


-65.7^1^:^ 


7 r:+40.6 
' -^-80.1 


84.3ti| 


n 9c;q+0.017 


0.3141°:°?° 


n 907+0.025 

^•^"^'-0.007 




rn 




6n 


O 


K 




6', a'n 




20.6111 


41.2t^^/ 


r-i 4+9.8 




'^••5^-0.00 




(0 ± 10) 


(0 ± 0.05) 



From Table [H by comparing the value and the number of data used in each group, we 
are able to have a rough idea on the quality of the fit. In most cases, the values are com- 
patible or smaller than the numbers of data used, indicating reasonable fit to measurements 
in these groups. However, the xIa{b-^K7t),...}^ xlAiB-^nw),...} ^^'^ ^5(50) values are larger 
than the corresponding numbers of data used. We will discuss more on the sources causing 
these sizable later. 

We give the fitted parameters in Table [Tll Uncertainties are obtained by scanning the 
parameter space with < Xmm + 1- The parameters consist of those in factorization 
amplitude (in the upper table) and of FSI (in the lower table). Values given in parenthesis 
are not fitted ones. We take 6', a' = (0 ± 10)° and k', ^' = ± 0.05 for estimation. 



We note that: (i) Most of our fitted va^ 



ues for hadronic parameters in factorization 



amplitudes agree with those usually used in 2M, |4ll, |42l]. However, the fit seems to prefer a 



small value of F^=^(0), which is at the lower end of the allowed region given in Eq. (1281) . 
(ii) Although it helps improve the fit, the effect of 0^;^ is sub- leading. On the other hand, the 
effect of (pA is prominent. The fitted (pA is around —66°, which is close to —55° as used in 
the so-called S4 scenario in QCDF [24]. When turning off FSI phases, our results should be 
similar to those obtained in the S4 scenario, (iii) The fitted r ~ 21° and u ~ 41° are closer 
to r = 24.1° and u = 35.3° of the exchange-type solution [see Eq. ( l25l) ] than to r = —41.8° 
and u = —19.5° of the annihilation-type solution [see Eq. ( 12^ ]. The exchange-type solution 
is more favorable, (iv) In this work, residual FSI is taken as left-over FSI that complements 
the FSI in factorization amplitudes. In principle, there is a possible double counting in (j)A,H 
and residual FSI phases. However, in practice the residual FSI is dominated by the exchange 
rescattering, which provides important effects on rates and CP asymmetries as we shall see 
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later. These effects cannot be easily obtained by varying pa,h and in reasonable ranges. 
In fact, numerically the x^/d.o.f. will not be reduced by freezing either of these parameters. 
Hence, both parameters are numerically important. 



A. Rates in and B Decays 

In Table llllt we show the CP-average rates of PP decays. In the table, Fac, 

"FSI" and FSI denote factorization, partial FSI and full FSI results, respectively. The FSI 
results are obtained with the best-fitted parameters shown in Table [TTl The factorization 
results are obtained by using the same set of the best fitted parameters, but with the residual 
FSI phases {5'^'\ cr^'\ k^'^ and ^*^'-') set to zero, while the partially FSI results are obtained 
similarly, but only with the real FSI phases {5^'\ o"*^')) set to zero. Recall that in Eq. fl20l) the 

1/2 

Sres matrix can be factorized into two parts, one involving real FSI phases and the other 
involving imaginary phases. The "FSI" results only make use of the one involving imaginary 
phases and are sort of "half-way" from the factorization results to the full FSI ones. 

Uncertainties for factorization results are not given and can be found elsewhere (for ex- 



ample, in 2^). The first uncertainties in FSI results are obtained by scanning the parameter 
space with ^ Xmin + while keeping 5' = a' = and k' = ^' = 0. The second uncertain- 
ties in FSI results are from the variations of 5', a', and From the table, we see that 
the B Kn and B Kt]^'^ rates are quite sensitive to 5', a', k' and ^' . Hence, a larger 
variation of these parameters is not preferred by the data. 

As shown in Table llllt the residual FSI results agree with data. Before turning on 
the residual FSI, the factorization results are close to the S4 ones as expected. After the 
residual FSI is turned on, some rates are enhanced remarkably. In particular, 5° decays in 
the AS* = transitions receive large contributions from the residual FSI. In the following, 
we will focus on effects of the FSI on some interesting modes. 

Through the residual FSI, B^ tt+tt^ and ti^ti'^ rates ^ are reduced and enhanced roughly 
by factor 2, respectively, leading to a better agreement with data. Note that in the "FSI" 
case, the tt+tt^ rate is enhanced, while the tt^tt" rate is slightly reduced. Both of them are 



2 



For the factorization amplitudes, we use the central values of Gegenbauer coefficients for the pion wave 



function, Q!j(2GeV) — 0.2 ± 0.1, used in 40| and do not consider the case of using a larger Gegenbauer 



coefficient, which leads to a larger tt^tt^ rate. 



15 



TABLE III: Branching ratios of various B ^ PP modes in units of IQ-^. Fac, "FSI" and FSI 

:"ull FSI results, respectively. See the maintext for details. 



denote factorization, partial FSI anc 
Experimental results are taken from 



Mode 




Exp 


Fac 


"FSI" 


FSI 




K---K+ 


19.4 ±0.6 


(16.0) 


(22.5) 


90 1+1.7+2.5 
^'-'•-'^-0.3-2.5 




Wtt^ 


9.8 ±0.6 


(7.2) 


(10.2) 


q 9+0.7+1.2 
^•^-0.2-1.2 




Wt] 


1.0 ±0.3 


(0.9) 


(1.7) 


-1 4+0.4+0.5 
-■■•^-0.1-0.4 






64.9 ±3.1 


(66.4) 


(62.3) 


Q+6-9 +9.2 
"J"J-y-10.6-8.1 


B- - 




23.1 ± 1.0 


(18.0) 


(26.1) 


99 c:+2.6+3.0 
^^•"J-1. 1-0.7 


B- - 




12.9 ±0.6 


(10.1) 


(14.3) 


T 9 /i+l. 5+1.6 
-■-^•^-0.2-1.6 


B- - 


► K-f] 


2.7 ±0.3 


(1.4) 


(2.5) 


9 1+0.6+0.6 
^•-^-0. 1-0.5 


B- - 


► K-7]' 


70.2 ±2.5 


(70.1) 


(65.0) 


7n S+6-6 +10.3 
'"-'•°-12. 3-9.2 


B- - 




o-oy_o.4o 


(5.18) 


(5.18) 


r: 1O+0.55+0.00 
^•J^°-0.38-0.00 


B~ - 




i.OU_Q 27 


(1.22) 


(1.77) 


1 4(^+0.35+0.15 
-"-•^"-0.04-0.13 


B~ - 


► TT~ri 


4.4 ± 0.4 


(4.10) 


(4.47) 


A 90 + 0.59+0.34 

^•^'^-0.23-0.37 


B- - 


► TT^rj' 


9 7+0.6 
^•'-0.5 


(3.09) 


(2.76) 


01 +0.19+0.65 
"^•"^-'^-0.51-0.54 






5.16 ±0.22 


(6.65) 


(7.56) 


r qn+1.92+0.39 

J->JU_Q 4g_Q 4Q 






1.55 ±0.35'^ 


(0.50) 


(0.36) 


1 r|4+0.12+0.10 
-■■■^^-0.55-0.08 




r]rj 


0.8±0.4(< 1.4) 


(0.21) 


(0.10) 


n 4ft+0.24+0.10 
'-'•^'-'„o.ll-0.08 




rjr]' 


0.5±0.4(< 1.2) 


(0.22) 


(0.24) 


n 00+0. 39+0.24 
U-oo_o40-0.21 




r]'r]' 


0.9l|]j(< 2.1) 


(0.16) 


(0.30) 


1 nR+1-16+0.36 
-'^•'^"-0.31-0.28 




K+K- 


•J-J^^-o.io 


(0.09) 


(0.05) 


n in+0. 35+0.10 
'J- '-'-0.02 -0.06 








(1.47) 


(1.56) 


1 in+0.46+0.12 
-•^•^^-0.12-0.11 






0.9±0.4(< 1.5) 


(0.26) 


(0.37) 


n 01 +0.05+0.06 
^^•"^^-0.01-0.06 






1.2 ±0.7^* 


(0.32) 


(0.22) 


n ^9+0.02+0.13 
^•^^-0.15-0.11 



"An S factor of 1.8 is included in the uncertainty. 
^An S factor of 1.7 is included in the uncertainty. 



pushed even further from the data. There are the real FSI phases {6, a) that will change 
these rates in the right direction. 
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FIG. 2: vr+vr-, 7r°7r° rates (left) and B° K+K', K^K° rates (right) versus 5 are plotted. 

The solid line is obtained by using all other parameters set to their best-fitted values, while the 
dashed line is obtained using the exchange-type solution for FSI parameters (see text). Bands are 
one-sigma ranges of experimental data. Theoretical uncertainties are not shown. Note that the 
fitted 6/Tr is around 0.3 (see Table . 

In Fig. [21 we show the tt+tt^ and ir^n^ rates versus S. The solid line is obtained by 

using all other parameters set to their best-fitted values, while the dashed line is obtained 
using the exchange-type solution for FSI parameters [see, Eq. fl25l) ] with r, z/ fixed, a = S 
and K = taken from the average of the central values of the fitted k, and ^. We see that 

—>■ TT+TT^ and Tr^Tr*^ rates are reduced and enhanced, respectively, as 6 is increasing. Both 
rates reach the measured ones at 5 ~ O.Stt. 

It is known that in order to have the vr^vr^ rate as large as observed, we need a sizable 
color-suppressed tree amplitude [43I]. In the residual FSI, a large color-suppressed tree 
contribution can be generated from the exchange rescattering. As shown in the upper part 
of Fig. m the color-allowed tree amplitude of the B^ tt+tt^ decay, a main FSI source in 
this sector, can produce a color-suppressed tree amplitude for the 5° — > tt^tt^ decay through 
the exchange rescattering. At the same time, the ir'^n' rate is reduced as it rescatters. 
We see that the exchange rescattering is responsible for the enhancement of tt^tt^ and the 
suppression of vr+vr". 

In Fig. [2], we show the B^ K^K~ and K^K^ rates. It is known that the K^K~ 
rate is sensitive to annihilation-type rescattering [2?! (corresponding to the and terms 
as depicted in Fig [Tt^c) and (d)). In the SU(3) case (solid line), for 5 < vr/2, the K^K~ 
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u,d u,d u,d u,d 



w- 

b K^'^^ s Q Q 




s s s s 



FIG. 3: Exchange rescattering in vr^vr^, B K tt^ and Bg rj^''^ri^'^ decays. 

constraint can be easily satisfied, while in the U(3) case (dashed line), the constraint on 5 
is even weaker. These features are understandable, since in both cases, the exchange-type 
rescattering, which cannot generate K^K~ final state by rescattering the tt+tt^ decay 

amplitude, is dominating. Note that the K^K^ rate is reduced through FSI, giving better 
agreement with data without violating the K^K~ bound. 

In summary, the residual FSI improves the agreement between theory and experiment for 
rates, in particularly, it resolves the discrepancy between data and theoretical expectations 
on n~^7r~ and Tr'^vr^ rates. 

B. Direct CP Violations in B^ and B^ Decays 

Results for direct CP asymmetries in B^, B^ PP decays are summarized in Table HVl 
In general, the residual FSI has a large impact on direct CP violations of many modes. In 
the following, we will focus on some interesting results. 

We first concentrate on the modes that lead to the Ktt puzzle in direct CP violation. 
We see that before the residual FSI is turned on (i.e. taking Sj-es = 1), we have A{B'^ — >• 
K^n^) ~ A{B^ K^n^) ~ —0.12 from the annihilation amplitude for 0a(— —66°). 
After turning on the residual FSI {Sres 7^ 1), the asymmetry ^(i?° K^tt^) changes from 
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TABLE IV: Same as Table IIIH except for the direct CP asymmetries A (in units of percent) in 
various B — > PP modes. 



Mode 




Exp 


Fac 


"FSI" 


FSI 






-9.8l:\i 


(-11.8) 


(-12.2) 


_q n+2-0+2.0 




Wtt^ 


-1 ± IS'' 


(3.3) 


(0.9) 


19 0+2.2+1.7 




K^r] 




(10.7) 


(2.1) 


9Q 7+8.0+3.3 

^O. / 1 Q 1 Q 




K^rj' 


4.8 ±5.1 


(0.2) 


(0.6) 


1 7+0.8+0.3 
1- ' —0.2—0.4 


B- - 




0.9 ±2.5 


(0.3) 


(0.2) 


n q+0.7+1.2 

>J.O_Q g_]^ J 


B- - 


► K-TT^ 


5.0 ±2.5 


(-11.8) 


(-10.3) 


A o+l. 4+1.9 
4-o_]^ 2-2.0 


B- - 


► K-T] 


-27 ±9 


(39.8) 


(28.2) 


97 q+8.6+10.8 

— o.U — D.o 


B- - 


► K-T]' 


1.6 ± 1.9 


(-2.6) 


(-2.6) 


q q + 1.0+0.5 

—0.5—0.5 


B- - 


► vr^vr'^ 


6±5 


(-0.06) 


(-0.06) 


n 01 —0 no 


B~ - 




12+^^ 

— 18 


(-3.5) 


(-1.8) 


1 9 0+9.1 +16.0 
— 12.8 — 17.8 


B^ - 


► TT~rj 


-16 ±7 


(19.7) 


(22.0) 


19 q+4.1+3.5 
— 2.9—3.2 


B- - 


► TT-rj' 


21 ± 15 


(22.8) 


(20.3) 


v:a 0+5.3 +1.7 
'-^ "— 10.6— 3.0 






38 ± 15'' 


(22.3) 


(21.1) 


1 c c;+10.2+4.6 
—4.3 —4.5 








(-51.5) 


(-45.8) 


AO q + 11.5+11.8 

"■"-*— 33.1 — 13.1 




r]rj 




(-11.7) 


(-77.6) 


crn 7+15.0 + 15.7 
— 12.4 — 16.3 




r]rj' 


— 


(-28.5) 


(-29.3) 


c- 7+9.5 +7.8 
"-••'-22.2-7.4 




ri'rj' 




(3.6) 


(18.7) 


9Q 7+26.2+8.3 
■^^•'-1.7 -6.6 




K+K- 




(0) 


(52.9) 


71 n+10.9+20.6 
' ^-^-41.4-15.6 






c:o+73 


(-9.0) 


(-19.9) 


q7 0+ 8.4+15.2 
•""•"-37.1-15.0 








(19.7) 


(19.1) 


7 9+11.5+0.4 
' -^-13.8-0.5 








(13.2) 


(12.4) 


99 7+ 7.7+1.0 
'-20.5-1.0 



"An S factor of 1.4 is included in the uncertainty. 
^An S factor of 2.4 is included in the uncertainty. 



~ —0.12 to ~ —0.09, while A{B~ — » K~ti^) changes from ~ —0.12 to ~ +0.05, reproducing 
the experimental results. In other words, the residual FSI modifies A{B^ K^ti^) and 
A{B~ — > K^Ti^) by an amount of ~ ±0.03 and ~ ±0.17, respectively. 
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-15 P ^ - 100 P ^ 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

5/7i: 5/ji 

FIG. 4: Same as Fig. [21 except that direct CP violations oiW K~Tr+ and B' K~7r° (left), 
and K~r] (right) versus 6 are plotted. Note that the best fitted iJ/vr is around 0.3. 

The residual FSI has a more prominent effect on A{B^ — > K^n^), and, hence, it is 
capable of lifting the degeneracy of A{B~ K^-k^) and A{B^ K^n^). As shown in 
Fig. m it only takes a small amount of 6 (~ 0.27r) to flip the sign of A{B^ K^tt^), but a 
large 6 (> O.Stt) would be needed to do the same thing on A{B^ — > K^n^). 

It is known that a sizable and complex color-suppressed tree amplitude in the B^ K^tt^ 
decay can solve the Ktt puzzle . As depicted in Fig. [3], a color-suppressed tree amplitude 
in the K~Tr^ mode can be generated from the exchange rescattering of B^ K^r]^'^ color- 
allowed tree amplitudes, which are known to be quite sizable [24,] . The rescattering leads to 
the desired large and complex color-suppressed amplitude in the K~it'^ mode. 

We note that in order to solve the Ktt direct CP puzzle, both (pA and 6 phases are needed. 
For example, a similar analysis using rescattering among naive factorization amplitudes that 
lack a large annihilation strong phase, was unable to remove the degeneracy of A{K^7t^) 
and A{K^7r^) In other words, rescattering from both in-elastic channels and PP final 
states contribute to A{K'n')s, reproducing the experimental results and resolving the Ktt 
direct CP violation puzzle without the need of introducing any new physics contribution. 

As noted in the previous section, the exchange rescattering is also responsible for the 
enhancement of the B^ tt^tt^ rate. In Fig. [5l we show a two-dimensional plot, exhibiting 
the correlation of the ratio B{B^ -k^-k^) / B{B^ tt+tt^) with the difference = 
A{B^ K^TT^) — A{B^ — s> K^TT^). The light shielded area is obtained by scanning over 
—71 < 6,a < TT and one-sigma ranges of r, z/, k and ^, while keeping all other hadronic 
parameters at their best-fitted values. The dark shielded area corresponds to the exchange- 
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0.3 

0.2 

0.1 
< 0.0 
-0.1 
-0.2 
-0.3 

FIG. 5: Correlation of the ratio B{W tt°tt^)/B{W vr+vr-) with the difference = 
A{K~TT^) — A{K~TT^). The light shielded area corresponds to the restricted SU(3) case, the 
dark shielded area corresponds to the exchange-type U(3) case and the solid line is the same as 
the previous one except keeping k = = —0.05. See the main text for more details. 

type U(3) case and is obtained by scanning over — tt < 6 = a < n, —0.35 < n = ^ < 0.35, 
while using r and u given in Eq. fl25|) and keeping all other hadronic parameters at their best- 
fitted values. The solid line is obtained in a similar manner except keeping k, = C, = —0.05, 
which is the average of the central values of the best fitted k, and ^. Note that in this 
case, only one FSI parameter 6 is varied. From the plot, we clearly see that the data 
can be easily reproduced and the exchange rescattering is responsible for generating sizable 
and complex color-suppressed tree amplitudes that account for the difference and the 
B{W 7r°7r°)/i3(5° n+n-) ratio at the same time. 

We now continue to discuss FSI effects on direct CP asymmetries. There are several 
interesting results and remarks: (i) Large effects of residual FSI on A for several other modes 
are obtained. Direct CP asymmetries in 5° K'^n^, K^r], 7r°7r° decays and in K^n^, 
K~~r], Tx~ri decays change signs in the presence of FSI. In Fig. HI we see that A{B^ K^rj) 
is quite sensitive to FSI. The solid line passes through the one-sigma range of data around 

5 ~ 0.37r. (ii) Recall that in Tabled we have X{^(_B-^E-,r) } ~ "^-^ ^^'^ 'x\A{B-^TnT) } ~ ^-^ 
from A{B- , K'n^, K'T], K^r]') and A{B- -n' -n^ , , -n' 7] , -n' r]') results, 

respectively. We see from Table HVl that the main contributions to these are from A{B~ — >■ 
K~ri') and A{B~ — * tt"//'), respectively, (iii) Note that the experimental uncertainty of 
A{B^ TT+TT^) is enlarged by a PDG ^-factor originated from two different measurements: 
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TABLE V: Same as Table IIIIl except for the branching ratios (upper table) in the unit of 10 ^ 
and direct CP asymmetries (lower table) in the unit of percent for various Bs — > PP modes. 



Experimental results are from 2 



38|. 



Mode 



Exp 



Fac 



"FSI" 



FSI 



13(Bs° - 




5.00 ± 1.25 


(4.72) 


(6.08) 


, _Li c;7_i_n 9n 
^•°-'^-0.39-0.22 


B{Bs^ - 


K^tt") 




(0.68) 


(0.59) 


1 1 n+n 24+n 05 

-"^•^"^-0.33-0.04 


B(Bs^ - 






(0.28) 


(0.21) 


^ — .n._i_n 10-1-0 f\A 
*^-^^-0.16-0.04 


B(Bs° - 


Wri') 




(2.33) 


(2.11) 


, , _j_n 1 A ifi A 9 

^•^^-0.44-0.36 


B(Bs° - 
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FIG. 6: Same as Fig. [2] except that Bs° i^"7r+ and K+K' rates (top), rj^'^r]^''^ rates 

(bottom) versus 5 (left) or cr (right) (with all other parameters fixed at the best-fitted values) are 
plotted here. Theoretical uncertainties are not shown. Note that the best fitted values for these 
FSI phases are 6/-K ~ 0.3 and cr/vr ~ 0.5. 



0.25 ± 0.08 ± 0.02 and 0.55 ± 0.08 ± 0.05 from BaBar 



37| and Belle [4J], respectively. Our 



fitted result of A{B° tt+tt" 



'15.5 



+10.2+4.6>| 
-4.3 -4.5J 



prefers the BaBar data, (iv) The direct 



CP violation of 71^71^ flips sign, resulting a large and positive A^h^ti^). (v) The direct 

CP violation of tt^tt^ is very small and does not receive any contribution from the 

residual rescattering, since it can only rescatter into itself [see Eq. (I7j)]. The smallness of 
A{B~ TT^vr") is consistent with a requirement followed from the CPT theorem 45|]. The 
A{B^ 7r^7r°) measurement remains as a clean way to search for new physics effects 22|. 



C. Rates and Direct CP asymmetries in Decays 



We now turn to B^ decays. In Table |Vl we show the CP-averaged rates and direct 
CP violations of B^ —>■ PP decays. The results are then compared with data. From the 
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table, we see that: (i) The B{Bs K~tt~^) rate agrees well with data. From Fig. [6l we 
note that the result is in agreement with data for < 5 < 7r/2 and any value of a. (ii) 
The B{Bs K^K^) rate plotted in Fig. [6] versus 5 and a agrees with data, (iii) The 
B{Bs TT+TT^) data can be reproduced, but the result has a large uncertainty, (iv) The 
A{Bs — * K'n^) data can be reproduced, but the fitted value is close to the lower end of 
the data. 

We expect the residual FSI to have sizable contributions to various Bs PP decay rates. 
For example, from Fig. [3] we see that the Bs — > rj^'^r]^'^ decays also receive contributions from 
the exchange rescattering. Plots of Bg — > rf^'^rj^'^ rates versus 5 and a are shown in Fig. [61 
The Bg 7]' 7]' rate is quite sensitive to the FSI phase a. As shown in Table IVl the 5° rj'ri' 
branching ratio is enhanced by a factor of 1.2 and reaches 1.0 x 10~^, which can be checked 
in the near future. The Bg K^tt^ and K^rj modes are also quite sensitive to the residual 
rescattering and their rates are enhanced by factors of 1.5 to 2, respectively. 

Similar to B^^d cases, the residual FSI also has large impacts on many A{Bs PP). 
From Table |Vl we see that signs of A{Bs A'°7r°) and A{Bs K^rj') are flipped. Note 
that direct CP asymmetries in Bg K^tt^, K^t], tt^i] and tt'^t]' decays are close to or greater 
than 50%. On the contrary, direct CP asymmetries in penguin dominated 6 — > s transition 
modes, such as Bg K^K^, rji]' and rj'i]' decays, are predicted to be quite small. It should 
be noted that these results may be subject to some small SU(3) breaking effects. 

There are increasing interests in the Bg sector Q, expected that more data from 

CDF and other detectors should be available soon. Predictions on Bg decay rates and direct 
CP violations given here can be tested in the near future. 



D. Time-dependent CP violations in and Decays 

Results on time-dependent CP-asymmetries S are given in Table IVI[ The sources of the 
first two uncertainties are the same as those in previous tables, while the last uncertainty 
comes from the variation of 7/03 = (67.6t4;5)° 129.]- We fit to data on mixing induced 
CP asymmetries. Note that for the B^ K^K^ mode, the mixing induced CP asymmetry 
obtained by BaBar {-l.2^tlf^tl^]^ [46]) and Belle (-0.38l[J:??±0.08 |47|) are quite different, 
and the central value of the former exceeds the physical range. Consequently, for this mode, 
only the Belle result is used in our fit. 
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TABLE VI: Results on the time-dependent CP asymmetry S of various B^^s PP modes. The 
first two uncertainties are same as those in previous tables, while the last uncertainty comes from 
the variation of 7/^3. 
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Time-dependent CP-asymmetries S of most decay modes, except —>■ vr'^vr^, rji], 
K+K~ and K^K° decays, do not receive large contributions from the residual FSI. Likewise, 
Sf in most of Bs modes are not sensitive to FSI effects, except those in B^ n^r]^'\ Kstc^ 
and Ksf]^'^ decays. 

For B^ decays, we define AS = sin 2/5efr — sin 2/3cc_ft:, where sin2/5efr = —rifS{f) with rjj the 
CP eigenvalue of the state /. Comparing with the recent value of sin 2j3ccK = 0.671±0.024 2I 
as measured in B^ K + charmonium modes, we obtain: 

AS{Ksn') = 0.107^^46, ^S{Ksv) = 0.098^°:°^^, AS{Ksv') = 0.011+°:°^6. (30) 

Note that the uncertainty in AS{Ksi]') are dominated by the one in the sin2j3ccK measure- 



found in 



HE 



ment. The AS{Ksr]'), being one of the promising tests of the SM {4], agrees with the one 
5|, while the ASlKs'n'^) given here is slightly larger. The main contribution to 



the xlsiW)} gi^^^ Table [T] is from S{B^ K^K^) 
We note that in the charming penguin approach 
CP violation puzzle AA{K'k) can also be resolved ISj and predictions on AS* are made. 



jy Ciuchini et al. [331], the Kn direct 



Ciuchini et al. obtained AS{Ksn^) = 0.024 ±0.059 and AS{Ksv') = -0.007 ±0.054 
Note that (i) their AS{Kstt^) overlaps with the one given in this work, (ii) the central value 
of their AS^Ksr]') is negative, but the associated uncertainties allow positive AS^Ksr]') as 
well. There is a considerable overlap between their AS{Ksri') and the one given here. 

For B'^ decays, the S contributed from B^-B^ mixing itself is around —0.036. Hence, for 
penguin dominated b ^ s transition, such as Bg K'^K^, rj^'^i]^'^ decays, we do not expect 
the corresponding \S\ to be much larger than (9(0.05). Indeed, the predicted \S\ as shown in 
Table IVT] for 5° —>■ rjr], ri'i]' and rjr]' decays are all below 0.06. In particular, given the large 
B^ rj^'^r]' rates, S{Bs r]'''^') are potentially good places to test the standard model. 
Given the recent interesting preliminary results in the Bg phase 
to search for S in these Bs charmless decays. 



23|, it will be very useful 



IV. CONCLUSION 

In this work, we study the FSI effects in all charmless -B«,d,s — ^ PP decay modes. We 
consider a FSI approach with both short- and long-distance contributions in which the 
former are from all in-elastic channels and are contained in factorization amplitudes, while 
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the latter are from residual rescattering among PP states. Flavor SU(3) symmetry is used to 
constrain the residual rescattering ^-matrix. We fit to all available data on the CP-averaged 
decay rates and CP asymmetries and make predictions on yet to be measured ones. Our 
main results are as follows: 

• Results are in agreement with data in the presence of FSI. 

• The fitted strong phase 0^ ^ —66° in annihilation amplitudes is close to the one used 
in the S4 scenario of the QCDF approach. 

• For B decays, the vr+vr" and n^-n-^ rates are suppressed and enhanced, respectively, by 
FSI. 

• The deviation (A^) between A{W K--n+) and A{B~ R- 7r°) can be understood 
in the FSI approach. Since A{K^Ti^) is more sensitive to the residual rescattering, the 
degeneracy of these two direct CP violations can be successfully lifted. However, both 
short and long distance strong phases are needed to give correct values for ^(i^7r)s. 

• It is interesting to note that the exchange rescattering is responsible for generating 
large and complex color-suppressed amplitudes, which are crucial in explaining the 
enhancements in the B{B^ ti^-k^) / B{B^ tt^tt^) ratio and the CP asymmetry 
difference A^. 

• The residual FSI has a large impact on direct CP asymmetries of many modes. 

• The direct CP violation of B~ ti~ti^ is very small and does not receive any con- 
tribution from the residual rescattering [see Eq. ([7])]. It remains as a clean mode to 
search for new physics phases. 

• The present data on Bg PP decay rates and direct CP violations can be successfully 
reproduced. 

• Several B^ decay rates are enhanced by FSI. In particular, the ri't]' branching ratio is 
predicted to reach 10~^ level, which can be checked experimentally. 

• Time-dependent CP asymmetry 5* in Bd,s decays are studied. The AS{B^ Ksi]') 
is very small (< 1%). This asymmetry remains as one of the cleanest measurements 
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to search for new physics phases. The fitted AS{B^ Ks7f^) is positive and cannot 
explain the present AS{B^ Kstc'^) data. 

• Most of the time- dependent CP asymmetries S of Bg to PP states with the strangeness 
S= +1 are expected to be small. The predicted \S\ for B^ rjr], rjr]' and rj'rj' decays 
are all below 0.06. These modes will be useful to test the SM. 
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APPENDIX A: MASTER FORMULA OF FSI 

Let ifw = J2q \0q denote the weak decay Hamiltonian, where Xq are VqbV*^ (or VqbV*^) 
and Oq are four-quark operators (including Wilson coefficients qs). From time reversal 
invariance of Oq, one has, 

{i;out\Oq\B)* = ((z;out|)*f/^f/TO*f/|f/T|5)* 
= {i;m\Oq\B) 

= in| k; out)(fc; out|Og|i?) 

k 

= ^5^,(fc;out|0,|E), (Al) 

k 

where Sik = (i; out| fc; in) is the strong interaction 5'-matrix element, and we have used 
?7T|out (in))* = I in (out)) to fix the phase convention. Eq. (1A1I) can be solved by (see, for 



example |48l ]) 



{z;ont\0,\B) = J2^}l'Af, (A2) 

k 

where y4|° is a real amplitude. To show that this is indeed a solution to Eq. flAll) . one 
needs to use Sik = Ski, which follows from the time reversal invariance of strong interactions 
and the phase convention we have adopted. The weak decay amplitude picks up strong 
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scattering phases [49|] and we have 



Af'^ ^ {t-oni\H^\B) = J2{r,ont\XMB) = Y.Sll\\Af) = (A3) 

q q,k k 

where we have defined = XqA'^^ free of any strong phase. The above equation is the 
master formula for FSI in Bu,d,s decays. 

APPENDIX B: CONSTRAINTS IN THE U(3) CASE 



In the U(3) case, one cannot have rescattering from both exchange and annihilation so 
that Te'^Vi'"'* = 0. This can be easily seen by inspecting 3 in the ti^-k^ — K^K^ — Ti^rjq — 
n^rjs basis, where 7]q = {uu + dd)/\/2 and ?7s = ss. From Eq. ([9]) and the requirement that 
iS^^ and T*^™'' preserve their forms as determined by U(3) symmetry, we should have 
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in the new basis. Under U(3) symmetry, it is evident that (5^^ 3)34,43 = for any m. Hence, 
from 



l'^res,3J34,43 — U + ^^^3 ^ ^3 ' ^3 J34,43 



which should also be zero, we must have 



r^Vi") = 



(B2) 



(B3) 



for any m in the U(3) case. 

Given the above constraint, we have two different solutions, which are, (a) annihilation 
type (ri™''' 7^ 0, ri™''' = 0) and (b) exchange type (re'"'' 7^ 0, ri™'' = 0). For the annihilation 
type solution, we have 



l+^(rr)+rM) 



ir 
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z(ri™)+r^)) 



20 



(B4) 
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which corresponds to taking 

c'27 = cis = cii = ^, ^8, di, T = --sm z/ = --sin (B5) 

in Eq. (fT5|) . For the exchange type solution, we have 

2m iJ _j_ g2m iSg ^ 









- ¥ 


' a 


— 

— 't 



which corresponds to setting 



(B6) 



027 = ^8 = ^1 = 5, ^8 = 5i, r = - sm — , z/ = - sm (B7) 

in Eq. ([T5|). 

In the above solutions, we explicitly see that U(3) symmetry imposes relations on the 
parameters of different SU(3) multiplets, and, consequently, reduces the number of indepen- 
dent parameters. It should be noted that mixing angles are fixed in both solutions. 
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